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ABSTRACT 

We report on a long (100 ks) XMM-Newton observation of the bright Seyfert 1 galaxy 
Arakelian 120. The source previously showed no signs of intrinsic reddening in its infrared- 
ultraviolet continuum and previous observations had shown no evidence for ionized absorp- 
tion in either the ultraviolet or X-ray bands. The new XMM-Newton RGS data place tight 
limits on the presence of an ionized X-ray absorber and confirm that the X-ray spectrum of 
Ark 120 is essentially unmodified by intervening matter Thus Ark 120 can be considered a 
'bare' Seyfert 1 nucleus. This observation therefore offers a clean view of the X-ray spectrum 
of a 'normal' Seyfert galaxy free from absorption effects. The spectrum shows a Doppler 
broadened iron emission line (FWHM ~ 3 x 10"* km s"^) and a smooth, continuous soft 
excess which appears to peak at an energy « 0.5 keV. This adds weight to the claim that 
genuine soft excesses (i.e. those due to a real steepening of the underlying continuum below 
^ 2 keV) are ubiquitous in Seyfert 1 spectra. However, the detailed shape of the excess could 
not be reproduced by any of the simple models tested (power-laws, blackbodies, Comptonised 
blackbodies, accretion disc reflection). This observation therefore demonstrates both the need 
to understand the soft excess (as a significant contributor to the luminosity of most Seyfert Is) 
and the inability of the existing, simple models to explain it. 

Key words: galaxies: active - galaxies: nuclei - galaxies: Seyfert - X-rays: galaxies 



1 INTRODUCTION 



Arakelian 120 (aka Mrk 1095) is a luminous Seyfert 1 galaxy at a 
redshift z — 0.0323. It was the subject of an early attempt at rever- 
beration mapping, and was important in demonstrating the compact 
size of the optical broad line region (Peterson et al. 1985; Peter- 
son & Gaskell 1991). More recent optical monitoring campaigns 
have yielded an estimate of the mass of the central black hole of 
~ 2 X 10* MoCWandel, Peterson & Malkan 1999). The bolo- 
metric luminosity for the nucleus is Lboi > 10^'' erg s~^ (Edel- 
son & Malkan 1996 estimated the total 0.1 — 100 nm luminosity 
to be ~ 8 X 10^* erg s^^). This would suggest it is radiating at 
L/L-Edd > 0.05, where L^dd is the Eddington luminosity for a 
2 X 10* Mq black hole. The nucleus is radio-quiet but shows a 
slight extension in its radio image (Condon et al. 1998; Ho 2002). 
Ward et al. (1987) used broad-band photometry to identify Ark 120 
as a 'bare' Seyfert nucleus, i.e. one free from significant reddening 
or contamination from the host galaxy. The host is a low-inclination 
spiral galaxy (Hubble type SO/a, inclination i « 26° ; Nordgren et 
al. 1995) . 



Ark 120 has been observed with most of the major X-ray ob- 
servatories. An EXOSAT observation showed Ark 120 to have a 
steep soft X-ray spectrum (Turner & Pounds 1989), as did a subse- 
quent i?05Ar observation (Brandt et al. 1993). Furthermore, these 
X-ray observations showed no indication of any 'warm absorption' 
features - i.e. discrete absorption features often found in the soft 
X-ray band and caused by absorption in photoionized gas along the 
line-of-sight to the nucleus. Warm absorption systems are common 
in Seyfert Is (Reynolds 1997; Crenshaw, Kraemer & George 2003). 
Observations in the ultraviolet (Crenshaw et al. 1999; Crenshaw & 
Kraemer 2001) showed Ark 120 to be one of the Seyfert 1 galax- 
ies that showed no intrinsic ultraviolet absorption. Thus Ark 120 is 
a rare example of a bright Seyfert 1 galaxy that is not significantly 
affected by any kind of complex absorption - its emission spectrum 
is that of a 'bare' Seyfert 1 nucleus. 

The paper presents the results of a long XMM-Newton obser- 
vation designed to characterise the intrinsic X-ray emission spec- 
trum of a Seyfert 1 galaxy. The rest of this paper is organised as 
follows. Section|2|discusses the observation details and data reduc- 
tion. Section |3| gives details of the X-ray variability observed in 
Ark 120 and section |4] describes a 'first look' at the X-ray spec- 
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trum. The spectrum is then examined in detail using the Reflection 
Grating Spectrometer (RGS) in section |5| This is followed by an 
analysis of the European Photon Imaging Camera (EPIC) data first 
over the 3—10 keV band (section|6j and then the over the full band- 
pass (section^. Section|8|briefly discusses an analysis of archival 
RXTE observations of Ark 120. Finally, the results are discussed in 
section|9| 



2 OBSERVATION AND DATA REDUCTION 

XMM-Newton (Jansen et al. 2001) observed Ark 120 during the pe- 
riod 2003 August 24 05:35:43 - 2003 August 25 12:44:33. XMM- 
Newton carries three co-aligned X-ray telescopes each with an 
EPIC CCD array as its focal-plane detector. There are two types 
of EPICs used: two MOS cameras (Turner et al. 2001a) and one 
pn camera (Striider et al. 2001). The two telescopes that focus onto 
the MOS cameras also focus onto RGS instruments (den Herder 
et al. 2001). During the Ark 120 observation all the EPIC cameras 
were operated in small window mode using the thin filter. The short 
frame-time possible in small window mode meant photon pile-up 
(which can distort the spectrum of bright sources; Ballet 1999) was 
negligible. 

The extraction of science products from the Observation 
Data Files (ODFs) followed standard procedures using the XMM- 
Newton Science Analysis System (SAS v5 . 4 . 1). The EPIC and 
RGS data were processed using the standard SAS processing chains 
to produce calibrated event lists. These removed events from the 
position of known defective pixels, corrected for Charge Transfer 
Inefficiency (CTI) and applied a gain calibration to produce a pho- 
ton energy for each event. Source data were extracted from circular 
regions of radius 45" from the processed images and background 
events were extracted from regions in the small window least ef- 
fected by source photons. These showed the background to be rela- 
tively low and stable throughout the observation. The total amount 
of 'good' exposure time selected was 78.2 ks from the pn, 108.4 ks 
from the MOS and 111.7 ks from the RGS. (The lower pn exposure 
is due to the lower 'live time' of the pn camera in small-window 
mode, ~ 71 per cent; Striider et al. 2001). The total number of 
source photons extracted was 2.2 x 10^ for the pn, 8.0 x 10^ for 
each MOS, and 8.6 x 10* and 9.9 x lO"* counts for the RGSl and 
RGS2, respectively. 

The spectral analysis was performed using XSPEC vll.2 
(Arnaud 1996) and, in order to allow effective use of the method of 
minimization for spectral fitting, the EPIC spectra were grouped 
such that each energy bin contains at least 50 counts. The RGS 
data were grouped to contain at least 20 counts per bin. The quoted 
errors on the derived model parameters correspond to a 90 per cent 
confidence level for one interesting parameter (i.e. a Ax^ = 2.71 
criterion), unless otherwise stated, and fit parameters are quoted for 
the rest frame of the source. Values of Ho = 70 km s~^ Mpc~^ 
and go = 0.5 are assumed throughout the paper. 



3 TIMING ANALYSIS 

Fig Q shows the broad band X-ray light curve obtained using the 
pn camera. This light curve has been background subtracted and 
corrected for exposure losses. The source shows only very mild 
variations in flux. Fitting the light curve with a constant gave an 
unacceptable fit (x^ = 929.8 for 109 degrees of freedom), reject- 
ing the constant hypothesis with very high significance (> 99.99 
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Figure 1. Top panel: EPIC pn light curve of Ai'k 120 (0.2 - 10 keV with 
1000 s bins). Note the false zero on the ordinate. The source clearly shows 
significant but low amplitude flux vaiiations. Bottom panel: softness ratio 
(0.2 — 0.5/2 — 10 keV) showing weak but significant spectral variations. 
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Figure 2. Rms spectrum of Ark 120. 



per cent) and demonstrating the observed variations are intrinsic to 
the source. The measured amplitude of variations, estimated using 
the Fvar statistic (see Vaughan et al. 2003), was only Fvar = 1.9 
per cent. 

The data were divided into four energy bands (0.2 — 0.5, 
0.5 — 1.0, 1 — 2 and 2 — 10 keV) in order to investigate the spectral 
dependence of the variations. The individual light curves show sim- 
ilar trends. The ratios of these light curves in the four bands were 
used to test of spectral variability. A constant was fitted to each 
of the ratios and the resulting value recorded. Spectral variabil- 
ity was clearly detected, with a constant hypothesis being rejected 
at > 90 per cent confidence in five of the six ratios, and rejected 
at > 99.99 per cent confidence for the ratio 0.2 - 0.5/2 - 10 
(X^ ~ 233.119/109 dof). The time series of this ratio is shown in 
Fig.Q The variability amplitude was measured in each of the four 
bands and used to define the rms spectrum (Fig.|2j. There was an 
obvious enhancement of variability amplitude at higher energies. 
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Figure 3. The 'fluxed' spectrum of Ark 120 after correcting for Galactic 
absorption. This was produced by calculating the ratio of the spectra of 
Ark 120 and 3C 273 and then normalising by the spectral model for 3C 273 
(in flux units). 
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Figure 4. The 'fluxed' X-ray spectrum of Ark 120 (see Figs.|5]and|5) and 

the non-simultaneous, multi-wavelength SED. 



4 A FIRST LOOK AT THE X-RAY SPECTRUM 

The X-ray spectrum in true flux units was estimated using the pro- 
cedure discussed in Vaughan & Fabian (2004). The counts to flux 
conversion was carried out using an XMM-Newton observation of 
3C 273 as a 'standard star.' 3C 273 is a bright source and has a rel- 
atively simple spectrum in the EPIC band (a hard power-law plus 
smooth soft excess modified by Galactic absorption; Page et al. 
2004a), and in particular does not contain any strong, sharp spec- 
tral features. The ratio of the spectra of Ark 120 and 3C 273 was 
calculated and then normalised by a spectral model for 3C 273 (see 
Vaughan & Fabian 2004 and also Page et al. 2004a). Fig|3|shows 
the resulting spectrum and clearly reveals the strongest spectral fea- 
tures in Ark 120, namely the smooth, featureless continuum and the 
iron line peaking at ^ 6.4 keV. 

This spectrum was corrected for Galactic absorption. The low 
Galactic latitude of Ark 120 (6 = —21.1°) means foreground 
Galactic absorption is high. The column density towards Ark 120 
has been estimated from 21 cm surveys of the surrounding sky to 
be N-a = 0.98 x 10^^ cm"^ (using the data of Stai'k et al. 1992) 
or Nh = 1.26 X 10^^ cm~^ (using the data of Dickey & Lock- 
man 1990). Either of these values can be consistent with the es- 
timate based on fits to the ROSAT soft X-ray spectrum (Brandt et 
al. 1993), although which one depends on the assumed form of the 
underlying spectrum. In the analysis below the Galactic absorp- 
tion was assumed to be given by the smaller of these two column 
density estimates (unless stated otherwise) and was modelled using 
the TBabs code of Wilms, Allen & McCray (2001). The relatively 
large correction for Ark 120 introduced spurious features into the 
fluxed spectrum at low energies. This is because the source spec- 
trum was convolved through the detector response but the absorp- 
tion model used for the correction was not. This lead to a spurious 
'inverted edge' feature appearing near the Galactic O I edge. Fur- 
thermore, the 3C 273 observation used the medium filter while the 
Ark 120 observation used the thin filter, this will systematically dis- 
tort the fluxed spectrum of Ark 120 at the lower energies. Due to 
these issues the fluxed spectrum is shown only above 0.8 keV. 

Figure |4] shows the X-ray spectrum compared to the over- 
all SED obtained from non-simultaneous observations. The multi- 
wavelength fluxes (VLA, NRAO, IRAS, optical) were obtained 



from the NASA/IPAC Extragalactic Database (NED). The ultra- 
violet HST spectrum was taken with the Faint Object Spectrograph 
(FOS) on-board HST and was previously published by Crenshaw 
et al. (1999). These data were obtained through the Multi-Mission 
Archive at Space Telescope (MAST) and de-reddened assuming 
E{B ~ V) = 0.128. The SED was produced using data obtained 
over a number of years by a number of observers and so should 
not be taken to represent the single-epoch SED of Ark 120. Nev- 
ertheless this SED does provide a strong indication of the general 
shape of the broad-band SED of Ark 120, which must peak in the 
far-ultraviolet region. Interestingly, the soft X-ray excess (see sec- 
tion0 does not appear steep enough to be the high energy tail of 
the blue-ultraviolet bump. 



5 RGS SPECTRAL ANALYSIS 

The combined RGSI-I-RGS2 spectrum is shown in true flux units 
(as a function of energy) in Fig|5| This was corrected for Galactic 
absorption using the TBabs model (which accounts only for pho- 
toelectric absorption) and assuming A'^h = 0.98 x 10^^ cm~^. Due 
to small inaccuracies in both the RGS response and the interstel- 
lar absorption model (see de Vries et al. 2003) this resulted in a 
spurious 'spike' in the corrected data at ~ 0.54 keV. Fig|6|shows 
a closer examination of the fluxed RGS spectrum as a function of 
wavelength. In this case no correction was made for the Galactic 
absorption and the only obvious features were those due to absorp- 
tion by interstellar O. Excepting this, the intrinsic soft X-ray spec- 
trum appears, to first order at least, as a smooth continuum with no 
obvious, strong, discrete features. 



5.0.7 Intrinsic warm absorption 

In order to test for traces of ionized (warm) absorption local to 
Ark 120 the RGS data were searched for resonance absorption 
lines of C, N, O and Ne. In particular, the (source frame) wave- 
lengths corresponding to the Is — 2p resonance transitions in He- 
and H-like ions of N, O and Ne as well as H-like C were examined. 
The relevant atomic data (wavelengths and oscillator strengths) 
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Figure 6. Close-up of the fluxed RGS spectrum plotted in wavelength units. The expected wavelengths of Galactic absorption features (the K-edge and 
principal resonance Hne of interstellar O I) are indicated as upwards arrows. Clearly the Galactic O I resonance line is detected. In addition, the wavelengths 
(shifted by the systemic velocity of Ark 120) of the principal resonance lines of highly ionized C, N, O and Ne, as well as neutral O, are shown with downwards 
arrows (all undetected). The slight 'dip' longwards of C VI coincides with a bad column in the RGS2 spectrum and is unlikely to be intrinsic to Ark 120. 
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Figure 5. Combined RGS spectrum shown in flux units after removing the 
Galactic absorption. (The dashed line shows the Galactic absorption profile 
scaled by a factor 2 X 10"'^.) Note that the emission line-like feature at 
^ 0.54 keV is largely an artifact of the absorption correction (compare 
with Fig|^. The nearby absorption line is real and due to the interstellar O I 
Is — 2p line. For clarity the error bars are not shown. 



were obtained from Verner, Vemer & Ferland (1996). The expected 
(source frame) wavelengths of these transitions are shown in Fig.|6| 

In order to measure the parameters of each line a region of the 
spectrum ~ 2 A either side of the expected transition was fitted 
with a power-law modified by Galactic absorption. A narrow (cr = 
0.1 eV) absorption line was then added to the data at the wavelength 
expected for the transition (after correcting for the redshift). In the 
case of all seven possible transitions no significant improvement 
in the fit was found upon adding an absorption line to the model. 
Upper limits were obtained on the equivalent widths of these lines 
as shown in Table Q These limits were derived after allowing for 
moderate red/blue shifts (up to 7000 km s~^) in the lines. Changing 
the wavelength range fitted did not substantially alter these results. 

These equivalent width limits place constraints on the possi- 
ble columns of ionized absorbing material as follows. Any weak 
absorption lines are likely to be unsaturated and thus on the lin- 
ear part of the curve-of-growth. For a given transition the measured 
equivalent width thus corresponds to an ionic column density given 
by the following (Spitzer 1978): 



X 



-^N,\f,j = 8.85 X lO^'^iVjA^- 



(1) 



where EW\ is the equivalent width and A is the wavelength 
(in cm) of the line, is the ionic column density and jij is the 
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Table 1. Equivalent width estimates for Is — 2p absorption lines as derived 
from the RGS data. The wavelengths are the expected rest-frame wave- 
lengths from Verner et al. (1996). " Galactic absoiption line (measured at 
z = 0.0). ' absorption line intrinsic to Ark 120 (measured at 2 = 0.0323). 



line 


A (A) 


EW{mk) 


Ni (cm^) 


C VI 


33.736 


< 64 


< 1.5 X lO^'' 


N VI 


28.787 


< 30 


< 6.1 X 10^5 


N VII 


24.781 


< 12 


< 5.3 X 10^5 


VII 


21.602 


< 28 


< 9.7 X 10^^ 


VIII 


18.969 


< 10 


< 7.5 X lO^s 


Ne IX 


13.447 


< 2 


< 1.7 X lO^s 


Ne X 


12.134 


< 12 


< 2.2 X 10^^ 


Ol° 


23.50 


60± 13 


6 - 10 X 10^'' 


Ol' 


23.50 


< 20 


< 1 X lOi'' 



oscillator strength. The derived limits on the ionic column densities 
are given in TableQ 



5.0.2 Interstellar absorption 

As mentioned above, the column density of Galactic interstellar gas 
is uncertain but estimated to lie in the range ~ 1.0 — 1.3 x 10^^ 
cm~^. The RGS spectrum clearly shows discrete structure due to 
absorption by local interstellar O I (see Fig.|6| note that the region 
of the spectrum around O I is sampled only by RGSl). The depth 
of the O I Is — 2p absorption line can provide an independent es- 
timate of the column density. The line is a simpler diagnostic than 
the O I K-edge which is broader and complicated by instrumental 
absorption (see de Vries et al. 2003). The equivalent width of the 
line {EW ~ 60 m A) corresponds to an absorbing column density 
log(iVo) ~ 17.6 - 18.0 (see Fig 3 of de Vries 2003). The main 
source of uncertainty in the column density estimate is the turbu- 
lent velocity width of the line; the line is deep enough that it does 
not lie on the linear part of the curve-of-growth. This column of 
neutral O corresponds to an equivalent hydrogen column density of 
Nn f» 0.8 - 2.0 X 10^^ cm"^ (assuming an O/H abundance of 
4.90 X 10^* as given by Wilms et al. 2000). The column density 
estimate is therefore consistent with the 21 cm measurements but 
does not improve upon them. 

The neutral absorption intrinsic to the host galaxy of Ark 120 
was estimated by measuring the O I Is — 2p line in the redshifted 
frame of the galaxy. The limit on its equivalent width is given in Ta- 
bleQ This corresponds to column densities of A'^o < 10^^ cm^'^ 
and TVh < 2 x lO^'' cm~^. The interstellar absorption intrinsic to 
Ark 120 is therefore negligible compared to the foreground absorp- 
tion. 



6 EPIC 3 -10 KEV SPECTRAL ANALYSIS 

The fluxed spectrum (section|4j demonstrated that the X-ray emis- 
sion from Ark 120 appears as a reasonably smooth continuum and 
the RGS showed no evidence for warm absorption (section|5j. This 
section describe an analysis of the EPIC spectrum, starting with an 
examination of the iron line. 




Observed energy (keV) 



Figure 7. EPIC pn and combined MOS spectrum of Ark 120. The data 
were fitted with a power-law over only the 3.5 — 5 keV 'continuum' region. 
The model was then extrapolated up in energy to reveal the excess due to 
the Fe emission. (NB: the MOSl and MOS2 data were fitted separately and 
combined only for plotting puiposes.) 



6.1 Iron Kq emission line 

The 3.5 — 10 keV spectrum was fitted, after excluding the 5 — 8 keV 
region which may contain Fe K features, with a simple model com- 
prising a power-law modified only by Galactic absorption. The pn 
and MOS spectra gave differing spectral slopes. The photon in- 
dex of the power-law fitted to the pn data was AF ~ 0.1 steeper 
than that fitted to the MOS data. The two MOS spectra were in 
good agreement with each other. Similar discrepancies have been 
seen in other observations (Molendi & Sembay 2003; Vaughan & 
Fabian 2004). Hence, in the following spectral analysis, the power- 
law slope was always allowed to differ between the pn and the MOS 
(but with the slope for the two MOS spectra tied together). 

The fluxed spectrum clearly showed Fe emission (Fig|3}, if 
this is produced by fluorescence in an optically thick medium then 
there may also be an associated Compton reflection continuum 
(Lightman & White 1988; George & Fabian 1991), which would 
show as a flattening of the spectrum at higher energies. Thus the 
spectral slope was estimated using only the 3.5 — 5.0 keV region 
- this is free from absorption effects and contains no significant 
contribution from either the line or reflected flux. The slope was 
measur ed to be r = 2.01 ± 0.06 (from the pn). FigQshows the 
residuals from this fit after the model was extrapolated to high en- 
ergies. The excess residuals due to Fe emission are obvious. There 
also appears to be a slight up turn in the continuum at energies 
above the line which could be due to reflection (see sections [6.21 
and[8j. 

The line could in principle be composed of lines at energies 
of 6.4 — 6.9 keV which could have profiles ranging from unre- 
solved (a < 70 eV) to very broad and asymmetric. A set of trial 
models were compared to the data to test the most plausible pos- 
sibilities for the line structure. The continuum was assumed to be 
a simple power-law between 3.5 — 10 keV (modified by Galactic 
absorption). The line was then fitted with a variety of alternative 
models including unresolved, narrow lines (modelled using narrow 
Gaussians of width cr = 10 eV), broad Gaussians and relativisti- 
cally broadened line profiles (the diskline model of Fabian et 
al. 1989). 

The best-fitting parameters and fit statistics are listed in Ta- 
ble H The seven different models tested were as follows: (1) one, 
(2) two and (3) three narrow Gaussians; (4) one broad Gaussian; (5) 
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Table 2. Best-fitting parameters for models of the iron line. Column (1) describes the line model components, either broad/narrow Gaussians or diskline 
profiles. Columns (2) and (3) give the line energies and equivalent widths. Column (4) gives the Gaussian width (cr) or emissivity index (q) in the case of 
diskline models. Columns (5) and (6) give the inner and outer radii and column (7) gives the incHnation angle for diskline components. Column (8) lists the 
fits statistics. ^ indicates the parameter was fixed. 
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(eV) 
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28 ± 5 


lO-^' 
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6.34 ±0.02 
6.76 ±0.02 


44 ±9 
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26 ±5 
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1580.0/1516 


4 


broad 


6.45 ±0.02 


130 ± 14 


177+11 








1596.0/1519 


5 


broad-l- 
nan'ow 


6.57 ±0.05 
6.40 ±0.01 


109 ± 16 
37 ± 7 


288t« 








1545.1/1517 


6 


diskline 


6.42 ±0.01 


156 ± 15 


-2.24 ±0.12 


""-39 


> 10^ 


> 48 


1566.3/1516 


7 


diskline-l- 
narrow 


6 56+0 02 
6.40 ±0.01 


42 ±6 


9 O-I+0.23 

lO-^' 


14434 


> 10^ 


> 73 


1541.8/1514 



line was described well in terms of either a Gaussian or a diskline 
profile. In the case of the broad Gaussian, the line width (a « 
290 eV) corresponds to a velocity width of FWHM « 3 x 10'' km 
s~^. This is far broader than the broad, permitted optical lines 
{FWHM{Uf3) « 5800 km s'^ Wandel et al. 1999) and so places 
the X-ray line emitting region within the optical Broad Line Region 
(BLR). Assuming the line originates from a Keplerian disk about a 
black hole the width measured from the broad Gaussian fit places 
the line emitting material within a few hundred Vg of the black hole 
(depending on the inclination). 

In the case of the diskline model the line profile was broad 
but reasonably symmetric - the large inner disk edge and high 
inclination angle mean Doppler shifts will dominate over gravita- 
tional redshift. The inner disc radius was constrained to lie beyond 
fin = 6rg, the innermost stable circular orbit (ISCO) expected for 
an accretion disc about a non-spinning (Schwarzschild) black hole. 
Fixing the inner radius to rin — 6r-g (in model 7) gave a signifi- 
cantly worse fit even with all the other line profile parameters left 
free (Ax^ = 14.0). That said, the x'^-space is complicated and 
smaller inner radii can be permitted provided the emissivity index 
(q where emissivity is J(r) oc r') is flatter and inclination angle is 
lower, as illustrated by Fig.|8| Nevertheless rin = 6r-g always lies 
outside the 99 per cent confidence region of the best-fitting model 
(model 7). 

The best-fitting energy of the line is rather unusual, at i5 ~ 
6.56 keV in both Gaussian and diskline models. A line energy of 
6.4 keV, corresponding to neutral iron, gives a substantially worse 
fit to the data (Ax^ = 21.5 compared to model 7, Table|2j. This 
best-fitting energy corresponds to the line emission being domi- 
nated by mildly ionized iron (Fe XX-XXII) meaning that resonant 
trapping and the Auger effect (see e.g. Ross, Fabian & Brandt 1996) 
should destroy the line. 
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Figure 8. Contours of Ax^ around the best-fitting diskline model (sec- 
tion l6. II as a function of disc emissivity index q and inner radius rin. The 
contours delineate the Ax^ = 2.3, 4.61, 9.21 regions. The best fitting so- 
lution (marked with a cross) has a large inner radius but the Ax^-space 
is complicated and lower inner radius solutions can be allowed provided 
the emissivity law is flatter. The other diskline pai'ameters (rest-frame line 
energy, outer radius and inclination angle) were free parameters. 



one broad and one narrow Gaussian; (6) a diskline; (7) a diskline 
and one narrow Gaussian. The data were not consistent with a 
model for the iron line comprising just one or two unresolved emis- 
sion lines. That said, the data did not strongly favour any one of the 
more complex models, models 3 — 7 all provided reasonable fits to 
the data (rejection probability < 0.95). However, the two models 
that include both a broad and a narrow line (models 5 and 7) gave 
superior fits. 

The best description for the iron line in Ark 120 is thus in 
terms of a narrow line (with equivalent width ~ 40 eV) and a 
stronger, broader line (with equivalent width ~ 100 eV). The broad 
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6.1.1 Additional iron line complexity 

The high inner radius of the diskline model can be alleviated 
if the intrinsically narrow line component is replaced by a second 
diskline. Fixing — 6 and rout = 10*rg for the pair of lines 
(with energies and emissivity indices free to vary independently) 
gave a good fit (x^ ~ 1545.8/1515 dof), comparable to the best 
fits given in table|2| The line energies and emissivity indices were 
as follows: Ex = 6.41 ± 0.01 keV, = -I-S+q.^ and E2 = 
6.80t°Q lo keV, /32 = -2.2~5';(^. The best-fitting inclination angle 
was i = 32 ± 5°. Thus it remains plausible that rin = 6 as long 
as there are two components to the Fe Ka line from the disc, each 
with flat emissivity laws. 

The only noticeable residuals that remained after the narrow 
and broad lines had been fitted (with either model 5 or 7), appeared 
as a slight excess at i5 ~ 7 keV and a slightly deficit at ~ 6 keV. 
The former of these could plausibly be accounted for by emission at 
6.9 keV by the Lya line of H-like Fe, or 7.06 keV emission the by 
the K/3 line of neutral Fe, or some combination of the two. Adding 
an additional narrow Gaussian to model 5 (Table|2j improved the 
fit by Ax^ = 9.0 with best fitting energy E = 7.02 ± 0.05 keV 
and an equivalent width EW = 15 ± 8 eV. The intrinsic width was 
poorly constrained when left free (a < 246 eV). The ratio of flux 
in the main Fe lines to the flux in the ~ 7 keV line was ~ 7.8 for 
the broad line and ~ 2.6 for the narrow line. The expected ratio of 
Ka/K/? flux is 150/17 ~ 8.8 for neutral iron. Thus the weak line 
at ~ 7 keV could be accounted for as the K/3 line accompanying 
the broader of the two Ka lines. However, some contribution from 
Fe Lya and/or a narrow, neutral K/3 line cannot be ruled out. 

The negative residual at i5 « 6 keV again seemed marginally 
significant. Including a narrow Gaussian absorption line improved 
the fit by a further Ax^ = 9.0 with best fitting energy E = 
6.03115:^4 keV and an equivalent width EW = -11 ± 6 eV. The 
identification of this feature is less obvious. As it is only a tenta- 
tive detection, rather weak (in equivalent width terms) and has no 
clear identification, a detailed analysis of this possible feature is not 
given in this paper. 



6.2 Reflection spectrum 

As noted above the X-ray continuum appears to be slightly curved 
compared to a power-law. This curvature is in the same sense as 
expected if there is a Compton reflection continuum component 
present in the data. Furthermore, an archival RXTE spectrum cov- 
ering the 3 — 20 keV range (see section [8} supports the presence 
of a reflection continuum. A simple hypothesis would be that the 
spectral curvature in the 3.5 — 10 keV EPIC spectrum is caused 
by a reflection continuum which originates in the same material as 
that responsible for the iron line. Assuming a neutral slab with so- 
lar abundances (George & Fabian 1991; Matt, Fabian & Reynolds 
1997) the equivalent width of the broad line (EW « 100 eV) cor- 
responds to a reflection strength of i? f» 0.5 — 1 (where R repre- 
sents the ratio of the normalizations of the observed and reflected 
power-law continua) depending on the viewing inclination. 

As a first test of this, a reflection continuum was included in 
the model using the pexrav code (Magdziarz & Zdziarski 1995). 
This component, which accounts for only the continuum produced 
by Compton reflection and photoelectric absorption (not any as- 
sociated fluorescent/recombination emission) was added to the ba- 
sic model of a power-law continuum plus two Gaussian iron lines 
(model 6 of Table|2j. The quality of the fit was only very slightly 
improved (x^ ~ 1542.2/1516 dof) and the relative reflection 



strength was constrained to lie in the range R — 0.02 — 0.89. This 
calculation was performed assuming the reflector was inclined at 
i — 30° with respect to the line of sight but almost identical con- 
straints were obtained assuming an inclination angle of i = 60°. 
Thus the XMM-NewtonfEPlC data are consistent with the presence 
of a moderate strength Compton reflection component (R < 1), 
which would match the strength of the iron line, but are not suf- 
ficient to give a significant detection (but see section[8j. The line 
parameters do not change substantially upon the inclusion of a re- 
flection component to the model. 

In order to test for the possible presence of a strongly gravita- 
tionally redshifted emission line this reflection plus double Gaus- 
sian model was refitted including an addition Laor emission pro- 
file (Laor 1991). The inner and outer radii of the disc were fixed at 
rin = 1.24 and rout = 20 to simulate only the very innermost re- 
gion of a disc extending close to the black hole. Assuming a neutral 
line and an emissivity power-law index of g = —3, the upper limit 
on the equivalent width of this feature was EW < 7, 21, 131 eV 
for inclination angles of i = 30° , 60° , 80° respectively. At high in- 
clinations, the line emission from this region of the disc is smeared 
over such a large energy range that, even in this relatively high qual- 
ity spectrum, it was difficult to constrain. It is clear however that 
unless the disc inclination angle is high i > 60° there can be no 
strong, gravitationally redshifted component to the line. 

The possibility of reflection was examined further using the 
reflection code of Ross & Fabian (1993) to model the line emis- 
sion and Compton reflection continuum self-consistently. The com- 
puted reflection spectrum was convolved with a diskline kernel 
to account for possible Doppler and gravitational broadening. The 
ionization of the reflector was kept fixed at its minimum allowed 
value (5 = 4nFx/nH = 10) and an additional narrow, neutral 
(6.4 keV) iron line was included using a (cr = 10 eV) Gaussian 
profile. This model therefore accounts for both the strength and 
profile of the broad iron line as well as the strength of the Compton 
reflection continuum. The best-fitting model provided an accept- 
able fit to the data (x^ = 1552.7/1518 dof) with parameters as 
follows: F = 1.96 ± 0.02, R = 0.49 ± 0.07, = 212ti33rg 
audi 43°. These results are consistent with both the diskline 
model discussed above (section l^TI and the RXTE fitting described 
in section[8| and show the broad iron line and weak reflection con- 
tinuum are consistent with a common origin. 



7 THE BROAD-BAND EPIC SPECTRUM 

The broad-band X-ray spectrum shown in Fig|3|is that of a smooth, 
curving continuum, closely resembling that seen in the X-ray spec- 
tra of some other Seyfert 1 galaxies including Ton SI 80 (Vaughan 
et al. 2002), PKS 0558-504 (O'Brien et al. 2001a) and Mrk 478 
(Marshall et al. 2003). In order to gain a better view of the shape 
of the broad-band spectrum, Fig|9|shows the data/model residuals 
when the simple model described in section l6m was extrapolated 
to lower energies (compare with Fig|3j. Clearly there is a strong, 
smooth soft excess which dominates at energies below 2 keV. A 
similar pattern of residuals is seen in all three EPIC spectra. This 
differs from the soft excesses seen previously in e.g. Ton SI 80 and 
Mrk 478 in that the excess appears to roll-over at energies below 
'-^ 0.5 keV (at least when compared to the underlying F ~ 2 hard 
X-ray model). This represents a 'characteristic energy' of the soft 
excess; the spectrum of the excess changes from steeper than F « 2 
above ^ 0.5 keV to a flatter spectrum at lower energies (see also 
Fig.|5j. A similar roll-over in the soft excess has previously been 
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Figure 9. Broad-band EPIC pn spectrum (top panel) and data/model resid- 
uals (bottom panel). The model spectrum (histogram) was fitted only above 
3.5 keV and subsequently extrapolated to lower energies. The model com- 
prises a power-law plus reflection continuum and two Gaussian emission 
lines (one narrow, one broad; see section l^j] . Galactic absorption was in- 
cluded in the model. 



seen in a few Seyfert 1 galaxies (e.g. Mrk 359, O'Brien et al. 2001b; 
NGC 4051, CoUinge et al. 2001). This cannot be accounted for in 
terms of excess neutral absorption, the curvature predicted by ab- 
sorption does not match the observed spectral form. The remainder 
of this section describes the results of fitting variety of continuum 
models to the data in an attempt to account for the soft excess spec- 
trum (see Vaughan et al. 2002 for more discussion of these soft 
excess models). 

For the purposes of this analysis the following assumptions 
were made. The Galactic absorption was modelled using TBabs 
with the colunm density allowed to vary in the range A^h = 
0.9 — 2.0 X 10^^ cm^^ and there was assumed to be no absorption 
intrinsic to Ark 120. The models tested all accounted for the hard 
X-ray continuum using a power-law (see section ISHl . In addition 
a reflection spectrum was included (using the pexrav code un- 
less stated otherwise) as was a two-component iron line (modelled 
using the double Gaussian model discussed in section lSHl . The rel- 
ative strength of the reflection (R) as well as the normalizations of 
the two iron lines were accepted as free parameters (but with the 
energies of the lines held fixed). This model acted as the 'baseline' 
model to which various continuum components were added in an 
attempt to model the soft excess. The differences in the spectral 
slopes between the pn and the MOS (see section l6m can cause 
problems when attempting to simultaneously fit broad-band con- 
tinuum models to both datasets. Therefore only the analysis per- 
formed on the pn data is explicitly described, but at each stage a 
'sanity check' was performed by fitting to the MOS data indepen- 
dently. 

The trial soft excess models were as follows: (1) an addi- 
tional, soft power-law, (2) a break in the hard power-law, (3) two 
breaks in the power-law, (4) a blackbody, (5) two blackbodies, 
(6) three blackbodies, (7) a bremsstrahlung thermal spectrum, (8) 
a thermal Comptonization spectrum (using the CompTT code of 
Titarchuk 1994), (9) a disc blackbody spectrum (using the di skpn 
code of Gierliiiski et al. 1999), (10) thermal plasma (using the 
mekal code of Kaastra 1992 and Liedahl, Osterheld & Goldstein 
1995); (11) ionized disc reflection and (12) ionized reflector plus 
bremsstrahlung (13) two separate ionized reflectors. For models 
11 — 13 the pexrav reflection model and broad Gaussian were 
replaced with the ionized reflection model of Ross & Fabian (1993) 



Table 3. Fit statistics for various broad-band spectral models. The models 
are listed as follows: pi = power-law: bknpo = broken power-law; 2bknpo 
= doubly broken power-law; bb = blackbody; bremss = bremsstrahlung: 
CompTT = thermal Comptonization; diskpn = disc blackbody; mekal = 
thermal plasma; rf93 = ionized reflection (from Ross & Fabian 1993). The 
models are described in more detail in the text. 





Model 


xVdof 


1 


pl+pexrav+pl 


3192.0/1522 


2 


bknpo+pexrav 


2921.5/1522 


3 


2bknpo+pexrav 


1863.7/1520 


4 


pl-(-pexrav-(-l xbb 


2483.6/1522 


5 


pl+pexrav+2 x bb 


1891.4/1520 


6 


pl+pexrav+3 X bb 


1869.5/1518 


7 


pl+pexrav+bremss 


2227.7/1522 


8 


pl+pexrav+CompTT 


1944.2/1520 


9 


pl+pexrav -(-diskpn 


2199.0/1530 


10 


pl+pexrav+mekal 


4291.2/1530 


11 


pl+rf93 


2656.5/1522 


12 


pl-(-2xif93 


2130.4/1520 


13 


pl-|-rf93+bremss 


1898.2/1521 


14 


pl-|-rf93+pl 


2073.2/1521 



including convolution by a relativistic kernel (using the Lao r code; 
Laor 1991) to model emission from an ionized accretion disc. Ta- 
ble|3|summarises the results of fitting these different models to the 
broad-band EPIC pn spectrum. In no case was the fit acceptable^ 
(-Proj > 99.9 per cent). 

7.1 Power-law continua 

The two-slope power-law models 1 (double power-law) and 2 
(singly broken power-law) were unable to account for the roll-over 
in the soft excess and produced very large data/model residuals in 
the soft X-ray band. The power-law models were able to fit down 
down to ~ 0.6 keV, the point at which the spectrum starts to roll- 
over, but below this energy the source spectrum flattens substan- 
tially. This is the case even after allowing for additional absorption; 
including extra absorption in the model did not alter the fit and 
could not account for the roll-over in the soft X-ray spectrum. Thus 
the observed roll-over in the soft excess must be intrinsic to the 
source and not an effect of line-of-sight absorption. However, the 
addition of a much flatter power-law slope at the lowest energies 
greatly improved the fit (model 3). In this model the hard contin- 
uum has a slope of Fhi — 2.16 down to Em — 1.6 keV where is 
steepened to Fmod ~ 2.46, accounting for the upturn in the spec- 
trum. Below Eio = 0.45 keV the spectrum becomes very hard, 
with Flo = —0.46, corresponding to the low energy roll-over in 
the soft excess. 

7.2 Thermal continua 

The models based around thermal continuum emission (mul- 
tiple blackbodies, Comptonised blackbody, disc blackbody or 
bremsstrahlung) provided comparably good fits and managed to ac- 
count for some degree of the roll-over. Unsurprisingly, given the 
lack of emission lines in the RGS spectrum, the mekal model 
gave a very poor fit using solar abundances. When allowed to 

^ For this reason the Ax^ = 2.7 intervals are not very meaningful as 
confidence intervals and so are not quoted. 
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be free the best-fitting abundance was zero, which essentially re- 
produced a thermal bremsstrahlung continuum. Model 6, shown 
in Fig 1101 provided the best fit of the thermal models, using 
three blackbodies to model the soft excess (with temperatures of 
kT = 102, 234, 673 eV). The broad-band (0.3 - 10 keV) unab- 
sorbed luminosity (predicted by this three-blackbody model) was 
Lx ~ 2.8 X 10'"' ergs-^ 

7.3 Ionized reflection features 

The spectrum was also fitted with ionized reflection models includ- 
ing relativistic smearing. For all the reflection-based models the 
disc emission was very centrally concentrated (e.g. rin < 1.4rg for 
model 13 and rin 2.2rg for model 14). In fact, the best-fitting 
emissivity index was often unreasonably steep (q < ~8) and so the 
fits were performed with q — —6 kept as a fixed parameter. Fur- 
thermore, the steep emissivity meant the fitting was insensitive to 
the Tout parameter, this was therefore kept fixed at rout ~ 400rg. 

It was not possible to model the entire soft excess in terms 
of emission from an ionized accretion disc (models 11 and 12). A 
single ionized reflector (model 11) was unable to fit the soft and 
hard bands simultaneously. The spectrum above 3 keV was well 
fitted with only weak ionization (section \6l2\ . in which case the 
reflector does not produce enough soft X-ray emission to account 
for any of the soft excess. Allowing for an additional, independent 
ionized reflector (with much higher ionization; ^ ^ 10'^) still gave 
a poor fit to the data (model 12). The primary reason for this failure 
is that the highly ionized reflector, while producing a strong soft 
excess, also produced very strong, broad spectral features in the Fe 
K-band. Thus the bulk of the soft excess cannot be explained in 
terms of ionized reflection (see also discussion in Vaughan et al. 
2002). 

Ogle et al. (2004) recently proposed that weak structures in 
the soft excess of the bright Seyfert 1 galaxy NGC 4051 can be 
explained by O VIII recombination emission (both the lines and 
radiative recombination continuum, RRC) from a relativistic ac- 
cretion disc. Models 13 and 14 work on a similar principle to the 
model of Ogle et al. (2004), in that the bulk of the soft excess is ac- 
counted for using a continuum component (thermal bremsstrahlung 
or power-law) and the relatively weak residual structures are fitted 
using emission from an ionized, relativistic accretion disc). These 
models gave quite competitive fits to the data, superior to the ion- 
ized reflection models without an additional soft X-ray continuum 
component. In each case the disc was only mildly ionized (^ < 15) 
and so the fits were performed with the ionization parameter kept 
fixed at 5 = 10. At this low level of ionization the reflector con- 
tributed modest equivalent width, broad emission features to the 
soft X-ray spectrum (mainly from recombination to 0®+ and O^^). 
For example, the equivalent widths of the emissions lines predicted 
by model 14 were ^ 36 eV and ^ 20 eV for O VII Hea and O VIII 
Lya, respectively. 

These models did provide reasonable fits to the data (com- 
pared to the other models) and so it is at least plausible that some 
spectral structure in the soft excess of Ark 120 could be accounted 
for by ionized disc emission. Model 13 accounts for the majority 
of the soft excess using a kT ~ 0.36 keV bremsstrahlung con- 
tinuum whilst model 14 used a F ~ 3.6 power-law. Notably the 
bremsstrahlung continuum provided a better fit to the soft X-ray 
spectrum than the power-law (see Table |3}. However, a problem 
with both these models is that the width of the observed Fe Ka line 
in over-predicted (see section l6m by the strong gravitational red- 
shift on the reflection spectrum (implied by the small inner radius). 
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Figure 10. Broad-band EPIC pn spectrum (top panel) fitted with model 
(6) and data/model residuafs (bottom panel). The model comprised a hard 
power-law plus reflection (pexr av) and iron line emission. The soft excess 
was accounted for with three blackbodies spread over a range of tempera- 
tures. The residual bump around 0.55 keV coincides with the instrumental 
and Galactic oxygen edges (a similar pattern of residuals was seen in the 
spectrum of the quasar 3c 273; Vaughan & Fabian 2004). 

It should be noted that the Ross & Fabian (1993) code was 
chosen to model the soft X-ray emission expected from an ion- 
ized disc, instead of fitting the O VIII emission independently as 
was performed by Ogle et al. (2004). The model used by Ogle et al. 
(2004) neglected a number of effects that must be important in a re- 
alistic disc. Firstly it did not account for emission lines from other 
species (e.g. O VII) nor from other ions (e.g. C, N, Fe) which should 
also contribute to the soft X-ray spectrum (e.g. Ballantyne, Ross & 
Fabian 2002). Secondly it did not include the associated Compton 
reflection continuum. Thirdly, the O VIII Lyman series line ratios 
were calculated in the optically thin limit, whereas a realistic accre- 
tion disc is expected to be optically thick to the lines. This means 
that, before they can escape the disc, photons from higher-order 
Lyman series lines (e.g. Ly/3 etc.) will be degraded into lower en- 
ergy photons from other series and Lya photons. Thus the higher- 
order Lyman lines should be strongly suppressed, compared to the 
optically thin case, and the emission dominated by the Lya line. 
Therefore, a realistic accretion disc spectrum would be expected to 
produce negligible emission in the higher order Lyman lines. These 
effects are all included in the Ross & Fabian (1993) model. (See 
also Dumont et al. 2003 for more details of line transfer effects in 
ionized discs.) 

7.4 Residual problems 

Even the best of the models discussed above gave a statistically 
unacceptable fit to the data. The most prominent fit residuals were 
confined to three regions: 1.7 — 2.5 keV, 0.5 — 0.7 keV and <. 
0.35 keV. The first set of residuals were almost certainly due to 
instrumental Si and Au features (see e.g. Vaughan & Fabian 2004). 
The 0.5 — 0.7 keV residuals could plausibly due due to low level, 
broad emission structures (see above), but they coincide with the 
instrumental and Galactic oxygen edges (at 0.543 keV). The RGS 
data showed that there were no strong, narrow emission/absorption 
lines at ~ 0.55 - 0.60 keV (20.7 - 22.5 A). The MOS spectral 
residuals showed noticeably less structure than the pn residuals in 
this region. Thus, given the present calibration uncertainties, it is 
not yet possible to unambiguously assess the shape and significance 
of these residuals. Apart from these two specifics, the remaining 
residuals were smaller than <C 5 per cent and so the best of the 
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Figure 11. Top panel: RXTE light curve of Ark 120 covering the period 
1998 Feb 24 to 1999 Mar 14. The arrow indicates the period of intensive 
monitoring (1998 Dec 15-18). Bottom panel: light curve obtained from the 
intensive monitoring. 
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Figure 12. RXTE spectrum of Ark 120 taken from the period of inten- 
sive monitoring (1998 Dec 15-18). The three panels show the data/model 
residuals using the following three models: a simple power-law (top), and 
power-law plus Gaussian (middle) and a power-law plus Gaussian and re- 
flection (bottom). All three models included Galactic absorption. 



above models could be considered reasonable fits to the data (even 
though the over-all values were poor). The slight down-turn in 
the residuals below 0.35 keV (e.g. Fig. llO> may indicate that even 
for the best-fitting models the spectral curvature was not enough to 
explain the spectral roll-over. 



8 RXTEJiKTA. 

Ark 120 has been monitored by RXTE for several years. Fig llll 
shows a section of the RXTE light curve spanning approximately 
one year which clearly reveals the source is variable but over 
timescales longer than sampled by XMM-Newton. Marshall et al. 
(2004) discuss the RXTE monitoring in more detail. 

Over the period 1998 Dec 15-18 RXTE observed Ark 120 
rather more intensively, resulting in 67.6 ks of useful data from 
3 PCUs. These data were reduced using standard methods (see 
e.g. Markowitz, Edelson & Vaughan 2003) and a time averaged 
3 — 20 keV spectrum was extracted. The spectrum clearly reveals 
the iron line and also suggests the presence of a reflected contin- 
uum. 

Three simple models were compared to the data as detailed in 
Table0 These were (1) a power-law, (2) a power-law plus Gaus- 
sian and (3) a power-law plus Gaussian and reflection continuum. 
The reflection continuum was modelled using the pexrav code. 
Galactic absorption was included in all three models. The iron line 
was modelled using a single broad Gaussian, which was adequate 
given the low spectral resolution of the RXTE PCA. Fig |12l shows 
the residuals from each model. The simple power-law provided an 
unacceptable fit. The fit became acceptable once a Gaussian was 
added, but was significantly improved by the addition of a reflection 
continuum. Although these RXTE data were taken almost five years 
prior to the XMM-Newton observation, the RXTE spectrum lends 
weight to the XMM-Newton results by confirming the presence of 
an iron line and supporting the existence of a weak (R ~ 0.5) re- 
flection continuum. 

Very similar results were obtained from an analysis of the time 
average spectrum of the long timescale monitoring (excluding the 
1998 Dec 15-18 intensive monitoring data). 



Table 4. Best-fitting parameters for models of the RXTE spectrum. 



Model 


Parameter 


xVdof 


power-law 


r = 1.91 


631.4/44 


power-law + 


r = 1.93 ±0.01 


51.9/41 


Gaussian 


E = 6.47 ± 0.04 keV 






cr = 260 ± 90 eV 






EW = 217±20eV 




power-law + 


r = 2.01 ±0.03 


24.0/40 


Gaussian -I- 


E = 6.50 ± 0.05 keV 






o- = 328 ± 90 eV 






EW = 224 ± 22 eV 




reflection 


K = 0.50 ±0.19 





9 DISCUSSION 

9.1 The X-ray properties of Ark 120 

This paper describes the results of a long XMM-Newton obser- 
vation of the luminous Seyfert 1 galaxy Arakelian 120. The X- 
ray emission from the source was only weakly variable on the 
short timescales probed by XMM-Newton (section |3j but archival 
RXTE data revealed 'normal' Seyfert 1 variability traits on longer 
timescales (section |8l cf. Edelson & Nandra 1999; Uttley et al. 
2002; Markowitz et al. 2003b). The X-ray spectrum showed a no- 
table absence of warm absorption features (section|5}- The lack of 
absorption means that these data represent the 'bare' X-ray emis- 
sion spectrum of a fairly typical Seyfert 1. 

The spectrum above ~ 3 keV can be explained using a fairly 
conventional spectral model (sections 16.11 and I6.2t comprising a 
power-law continuum (F « 2) plus Doppler broadened emission 
from the surface of a weakly ionized reflector (with relative reflec- 
tion strength R ~ 0.5). The RXTE data support this model (sec- 
tion [SJ. The spectrum below ~ 3 keV becomes dominated by a 
steep, smooth, broad soft excess component which appears to peak 
at ~ 0.5 keV. 



9.2 The iron line of Arlt 120 

The resuhs of the spectral fitting indicated an iron emission line 
that is the composite of a weak, narrow line originating in distant 
material and a stronger line showing significant Doppler broaden- 
ing. 

A narrow, neutral iron line appears to be ubiquitous in Seyfert 
1 galaxies (e.g. Yaqoob, George & Turner 2002; Page et al. 2004b; 
Yaqoob & Padmanabhan 2004) and may have an original in the op- 
tical broad-line region or the putative molecular torus (see discus- 
sion in Yaqoob et al. 2001; Reeves et al. 2004). The lack of neutral 
absorption intrinsic to Ark 120 (section|5) confined this material to 
lie out of the line of sight, which in turn implies the covering frac- 
tion of the line emitting material must be below unity. The relative 
strength of the narrow line (EW ~ 40 eV) further suggests the 
sky covering fraction of the material is small (/c <J 0.5; Nandra 
& George 1994) or the optical depth is small (r < 0.1; Leahy & 
Creighton 1993) or both. 

The broad component to the line has a velocity width 
FWHM ~ 3 X 10* km s~^, far broader than the broad optical 
lines (e.g. FWHM(Ul3) « 5800 km s"^; Wandel et al. 1999). 
However, there is no requirement for the line emitting region to ex- 
tend into the region of strong gravity about the black hole (< 20 rg) 
which would produce an asymmetric redward tail on the line pro- 
file (section lSm . The best-fitting parameters of the disc line model 
are slightly unusual; the inclination angle is high (i > 73°), the 
inner radius is rather large (rin ~ 140rg) and the rest-frame en- 
ergy is unusual (E ~ 6.56 keV). One plausible origin for this line 
is a weakly ionized accretion disc. The fact that the inner radius 
is greater than 6rg could mean the disc is truncated. However, the 
unusual energy of the line (corresponding to mildly ionized iron) 
might allow for a simpler alternative. If the disc survives down to 
the ISCO but rapidly becomes ionized with decreasing radius then, 
depending on the detailed ionization structure, the innermost re- 
gions may produce little observable line emission (see e.g. Ross, 
Fabian & Young 1999). 

However, as demonstrated by Fig. |8| there is sufficient lever- 
age in the fit to allow for a smaller inner radius (rin <; lOOrg), 
and also lower inclination angle, provided the emissivity law is 
quite flat {q ~ —2). As mentioned in section lOI a good fit can 
be obtained with rin ~ 6rg and i ~ 30° provided that there are 
two co-existing Fe Kq emission lines arising in the disc (at ener- 
gies of 6.4 and 6.8 keV) and no intrinsically narrow component. 
This may be feasible if the disc is clumpy or inhomogeneous (Bal- 
lantyne. Turner & Blaes 2004). The required flat emissivity can 
be produced if the disc is illuminated from a great height. Fig- 
ure 1 131 shows the effect of altering r/h in a 'lamppost' geometry 
[i.e. the disc was assumed to be flat and illuminated by a central 
point source raised by a height h, in which case the emissivity law 
is given by J(r) oc h/{r^ + h^)^^^]. This demonstrates that if the 
X-illumination is provided from a height h ~ lOOrg the emissivity 
will be quite flat out to similarly large radii. Illumination from a 
large height would require a non-standard geometry for the X-ray 
emitting region, one possibility is the 'aborted jet' model of Ghis- 
elUni, Harrdt & Matt (2004). 



9.3 Ark 120 as a 'bare' Seyfert 1 

These new X-ray observations showed Ark 120 to posses no evi- 
dence of an X-ray warm absorber and placed upper limits on the 
ionic column densities that are substantially lower than those of 
more typical, absorbed Seyfert Is. For example, NGC 3783 shows 
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Figure 13. The local power-law slope of the emissivity function ( J(r) 
r?) for a lamppost geometry. 

O VII and O VIII absorption with column densities two orders of 
magnitude higher (Kaspi et al. 2002). Kaastra et al. (2000) observed 
NGC 5548 and measured C VI, N VI, O VII and O VIII absorption 
lines with corresponding ionic column densities an order of mag- 
nitude higher than the limits obtained for Ark 120. A last counter 
example is IRAS 13349-1-2438 (Sako et al. 2001) which shows ab- 
sorption fines of C VI, N VI, O VII, O VIII, Ne IX and Ne X, with 
column densities an order of magnitude higher than Ark 120. 

These X-ray observations therefore confirm that Ark 120 is 
indeed a 'bare' Seyfert 1 nucleus, as suspected based on its broad 
band spectral energy distribution (SED; Edelson & Malkan 1986; 
Ward et al. 1987) and the lack of ultraviolet absorption seen in 
HST/FOS spectrum (Crenshaw et al. 1999). Other ultraviolet obser- 
vations have been made with the Goddard High-Resolution Spec- 
trograph (Penton, Stocke & Shull 2000) and with FUSE (Wakker 
et al. 2003). Using these data Wakker et al. (2003) showed that 
the line-of-sight through the Galaxy towards Ark 120 possesses a 
very low column of O VI. A preliminary examination of these data 
suggest there may be weak traces of ionized absorption intrinsic to 
Ark 120. 



9.4 Other bare Seyferts 

Shortly after the original discovery of the X-ray warm absorber 
(Halpern 1984), this phenomenon was found to be common in 
Seyfert 1 galaxies. Low resolution spectra from EXOSAT and 
Ginga suggested ~ 50 per cent of Seyfert Is showed evidence 
for warm absorption (Turner & Pounds 1989; Nandra & Pounds 
1994). Subsequent observations with the better spectral resolution 
afforded by ASCA suggested the incidence of warm absorbers was 
possibly higher (~ 50 — 70 per cent; Reynolds 1997; George et al. 
1998). Similar results were found for the ultraviolet absorber based 
on HST FOS spectra (Crenshaw et al. 1999). 

The increased sensitivity and resolution offered by XMM- 
Newton and Chandra have allowed for even more sensitive searches 
for absorption. These have confirmed that there remains a signifi- 
cant population of Seyfert Is that lack a strong X-ray warm ab- 
sorber. In addition to Ark 120, other Seyfert Is that show a distinct 
lack of warm absorption in both their HST FOS ultraviolet and 
XMM-NewtonI Chandra X-ray spectra include Mrk 478 (Marshall 
et al. 2003a), Mrk 335 (Gondoin et al. 2002), Fairall 9 (Gondoin 
et al. 2001) and Mrk 205 (Reeves et al. 2001). All four of these 



12 Vaughan et al. 



were in the sample of Crenshaw et al. (1999) and classified as hav- 
ing no intrinsic ultraviolet warm absorber. Other objects known to 
lack X-ray absorbers include PKS 0558-504 (O'Brien et al. 2001a), 
MCG-2-58-22 (Weaver et al. 1995) and Ton SISO^ (Turner et al. 
2001b; Vaughan et al. 2002). 

These objects are probably not completely without ionized ab- 
sorption systems. Kriss (2002) reported weak absorption by the 
O VI AA1032, 1038 A resonance doublet m both Ton S180 and 
Mrk 478 based on high resolution FUSE data (see also Turner et 
al. 2002 for the FUSE observation of Ton SI 80). Such absorption 
is too weak to have been detected in the Crenshaw et al. (1999) ul- 
traviolet survey. The corresponding X-ray warm absorbers in these 
objects may be present but so weak as to have a negligible effect on 
the available data. 

The reason for the lack of an X-ray warm absorber is unclear. 
The most obvious explanation is that overall colunm density is 
lower in these objects (by at least an order of magnitude) compared 
to more typical, absorbed Seyferts. However, it is also plausible that 
a similar column of ionized gas exists but is either too highly ion- 
ized to show significant spectral features or lies out of the line of 
sight (requiring a covering fraction less than unity). A detailed sur- 
vey comparing the X-ray and ultraviolet emission/absorption line 
spectra of a large sample of bright Seyferts may be able to an- 
swer this question (see discussions in Crenshaw et al. 2003). A 
more speculative solution, recently proposed by Gierlinski 8l Done 
(2004), is that objects that lack narrow absorption lines may be 
dominated by a deep, broad absorption trough produced by an ab- 
sorption system with such high dispersion velocity that no individ- 
ual lines can be resolved. This model would require the underlying 
power-law continuum to be rather steep, which then makes explain- 
ing the upturn above ~ 10 keV more difficult to explain. 

9.5 How common are soft X-ray excesses? 

These observations have clearly revealed a luminous soft excess in 
Ark 120, in the sense that an extrapolation of the hard (> 3 keV) 
continuum into the soft X-ray band revealed a very significant up- 
turn. Steep soft X-ray spectra were first seen in Se3^erts using 
HEAO-1 A2 data (Pravdo et al. 1981) and subsequently seen as 
excesses over the hard power-law by EXOSAT (AxndiViA et al. 1985) 
and Einstein (Bechtold et al. 1987). Early soft X-ray surveys (e.g. 
the £:X05Ar survey of Turner & Pounds 1989) suggested that > 50 
per cent of all Seyferts possessed a soft excess component. The 
/?05Ar survey of Walter & Fink (1993) also suggested a high inci- 
dence of soft excesses. These numbers were highly uncertain, how- 
ever, because many of the sample members were absorbed. Com- 
plex soft X-ray absorption can mask, or in some cases even mimic, 
the steeper soft X-ray spectrum indicative of a soft excess. 

The recent XMM-Newton and Chandra observations would 
seem to indicate that all Seyfert Is without strong X-ray warm ab- 
sorption do show a strong soft excess. However, many of the well- 
known examples (e.g. Ton SI 80, PKS 0558-504, Mrk 478, NGC 

2 There is some debate about Ton SI 80. R6zatiska et al. (2004) re- 
examined the Chandra LETGS observation and claim to have found several 
resonance absorption lines. The lines are very weak and in most cases of 
borderline significance (many have equivalent widths consistent with zero). 
Furthermore, the spectral analysis was complicated by the presence of con- 
taminant on the Chandra ACIS which affects the LETGS/ACIS calibration 
(see Marshall et al. 2003b). Thus the existence of an X-ray warm absorber 
in Ton SI 80 is still not well established and in any case must be extremely 
weak. 



4051 and Mrk 359) are narrow-line Seyfert Is (NLSls; Osterbrock 
& Pogge 1985). NLSls are a subclass of Seyfert Is defined by 
their narrow permitted optical lines (FWHM (H/?) < 2000 km 
s~^) but noted for their often exceptionally steep soft X-ray spec- 
tra (Boiler et al. 1996; Laor et al. 1997; Vaughan et al. 1999; 
Leighly 1999). Until recently it was possible that the soft excess 
appeared ubiquitous in unabsorbed Se3^erts only because many of 
the well-studied examples were NLSls (i.e. the sample of unab- 
sorbed Seyferts was biased towards the soft NLSls). 

Ark 120 is an interesting counter example, being both 
unabsorbed and a 'normal,' broad-line Seyfert 1 (BLSl) with 
FWHM{Yip) « 5800 km s"^ (Wandel et al. 1999). Yet it too 
shows a strong soft excess. Other notable BLSls that also lack ab- 
sorption include Mrk 205 (Reeves et al. 2001), Mrk 335 (Gondoin 
et al. 2002) and Fairall 9 (Gondoin et al. 2001), all of which showed 
soft excesses in their XMM-Newton observations. Although a com- 
plete and unbiased survey of the soft X-ray spectra of Seyfert Is 
has yet to be conducted, it does seem highly plausible that all unab- 
sorbed Seyfert Is (whether NLSl or BLSl) possess a soft excess. 
Asstuning that there is no fundamental difference in the underly- 
ing X-ray continuum spectra of absorbed and unabsorbed Seyfert 
Is then implies that soft X-ray excesses are ubiquitous to Seyfert 
Is. The mini-survey of XMM-Newton spectra by Pounds & Reeves 
(2002) showed that in all six Seyfert Is they studied, the soft X- 
ray spectrum (< 0.5 keV) was always > 50 per cent higher than 
an extrapolation of the hard X-ray power-law would predict. Fur- 
thermore, detailed studies of Seyfert Is with complex, strong X- 
ray warm absorbers often conclude that an additional soft X-ray 
emission component is required behind the warm absorber (e.g. 
CoUinge et al. 2001; Netzer et al. 2003; Blustin et al. 2003). 

This therefore underlines the need to understand the soft ex- 
cess as a common (perhaps ubiquitous) and significant contributor 
to the luminosity of Seyferts. In the case of Ark 120 the best-fitting 
models are a doubly-broken power-law with an anomalously flat 
slope below 0.5 keV, or multiple, soft blackbodies. The blackbody 
origin is difficult to explain as the temperatures are far too high to 
correspond to any standard accretion disc. The expected tempera- 
ture for the inner region (~ 6rg) of a geometrically thin, optically 
thick disc about a 2 x 10® M© black hole is feT ~ 11 eV if radiat- 
ing atL/LEdd = O.landfeT 20 eV if radiating atL/LEdd = 1. 
This emission should therefore not contribute to the observed X-ray 
spectrum. The temperatures of the best-fitting blackbody compo- 
nents were at at least an order of magnitude higher than this. In 
addition, the size of the emission region implied by these high tem- 
peratures is far too small (reB <C rg). The alternative models tested 
(reflection, disc blackbodies and bremsstrahlung) all produced the 
wrong spectral shape. At present there is no single model that can 
accoimt for the known properties of the soft excess (see also the 
discussion in Vaughan et al. 2002). 
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